Spins coupled to a Spin Bath: Prom Integrability to Chaos 
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Motivated by the hyperfine interaction of electron spins with surrounding nuclei, we investigate 
systems of central spins coupled to a bath of noninteracting spins in the framework of random 
matrix theory. With increasing number of central spins a transition from Poissonian statistics to the 
Gaussian orthogonal ensemble occurs which can be described by a generalized Brody distribution. 
These observations are unaltered upon applying an external magnetic field. In the transition region, 
the classical counterparts of the models studied have mixed phase space. 



Spins coupled to a bath of other spin degrees of free- 
dom occur in a variety of nanostructures including semi- 
conductor quantum dots carbon nanotube quan- 
tum dots 5j, phosphorus donors in silicon |6], nitrogen 
vacency centers in diamond , and molecular magnets 
|10| . A large portion of the presently very high both ex- 
perimental and theoretical interest in such systems is due 
to proposals to utilize such structures for quantum infor- 
mation processing [Tmi3j . Here the central spins play 
the role of the qubit whereas the surrounding bath spins 
act as an decohering environment. In the present let- 
ter we investigate very basic properties of such so-called 
central spin systems in terms of spectral statistics and 
random matrix theory 14 . 

The generic Hamiltonian is given by 

N 

n = Y.s^-J2A^h: (1) 

a—l i—1 

describing Nc central spins Sa coupled to A'' bath spins li, 
typically N ^ N^. Here we take all spins to be dimen- 
sionless quantum variables such that the coupling con- 
stants A^"^ have dimension of energy. A paradigmatic 
example is given by, say, a single spin of a conduction- 
band electron residing in a semiconductor quantum dot 
and being coupled via hyperfine contact interaction to 
the bath of surrounding nuclear spins. In a very typical 
material like gallium arsenide all nuclei have a spin of 
7 = 3/2 whereas in other systems like indium arsenide 
even spins of length / = 9/2 occur. In fact, this hyperfine 
interaction with surrounding nuclei has been identified 
to be the limiting factor regarding coherent dynamics of 
electron spin qubits [TH31 [H] . In the above example the 
hyperfine coupling constants A^'^^ are proportional to the 
square modulus of the electronic wave function at the lo- 
cation of the nucleus and can therefore vary widely in 
magnitude. For the purposes of our statistical analysis 
here we shall take an even more radical point of view 
and choose the a["'' at random. To be specific, we will 
choose the from a uniform distribution within the 
interval [0, 1] and normalize them afterwards according to 
J2i A^"^ — 1 for each central spin. The data to be pre- 
sented below is obtained by averaging over typically 500 



random realizations of coupling parameters. Note that 
the Hamiltonian matrix represented in the usual basis of 
tensor-product eigentstates of S^, /f is always real and 
symmetric. Therefore, the natural candidate for a ran- 
dom matrix description of such systems is the Gaussian 
orthogonal ensemble (GOE) [Til. 

In the important case of a single central spin, Ac = 1, 
the above model has the strong mathematical property 
of being integrable [inillT]. Moreover, this integrability 
is particularly robust as it is independent of the choice 
of the coupling parameters a[^^ and the length of the 
spins which can even be chosen individually 1 161 . In 
fact, the model ^ for a single central spin has been 
the basis of numerous theoretical studies on decoherence 
properties of c^uantum dot spin qubits; see, for example, 
Refs. [mils 21 , for reviews also p2H24] . It is an inter- 
esting question, both from a practical as well as from an 
abstract point of view, to what extend the results of these 
investigations are linked to the integrability of the under- 
lying idealized model. In particular, what changes may 
occur if the Hamiltonian deviates from the above simple 
case Nc — 1 by, e.g., involving more than one central 
spin? Previous investigations of decoherence properties, 
making strongly restrictive assumptions on the coupling 
constants, predicted a significant dependence on whether 
the number of central spins is even or odd [351 US] • In 
the following we will investigate Hamiltonians of the gen- 
eral type ([!]) within the framework of level statistics, i.e. 
generic spectral characteristics [13] ■ For other studies 
of interacting quantum many-body systems using this 
method see e.g. Rcfs. [27 29 . 

The spectra generated numerically from the Hamilto- 
nian ([T]) clearly have a nontrivial overall structure, i.e. 
the locally averaged density of states is not constant as 
a function of energy [l9l [22] . Therefore an unfolding 
of these spectra has to be performed which results in a 
transformation onto a new spectral variable s such that 
the mean level density is equal to unity [H]. We have 
compared several standard numerical unfolding proce- 
dures and made sure that they yield consistent results. 
Fig. [l] shows the probability distribution p{s) for the 
nearest-neighbor level spacing for a system of a single 
central spin 5*1 = 1/2 and 13 bath spins of length 1 — 1/2 
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for several subspaces of the total angular momentum 
J = Si + J2i where each multiplet is counted as a 
single enery level. The subspaces of highest J = 5, 6, 7 
have been discarded, and in the bottom right panel all 
probability distributions are joined. As to be expected 
for an integrable model, the level statistics follow a Pois- 
son distribution resulting in an exponential level spacing 
distribution p{s) = e~*. This is in contrast to the case of 
two central spins 5i = = 1/2 shown in Fig. [2) Here 
level repulsion takes clearly place, p(0) = 0, although the 
data considerably deviates from the Wigner surmise for 
the GOE im, p{s) = (7r/2)sexp(-(7r/4)s2). 

Obviously, our numerical studies are technically re- 
stricted to rather small system sizes, Nc + N < 14. This 
limitation, however, does not affect our results for the 
level statistics as demonstrated in Fig. [3] where we have 
plotted the same data as in the bottom right panel of 
Fig. [2] but for = 10 and iV = 11 bath spins. This 
insensitivity to the system size seen in the figure is a 
natural consequence of the unfolding of the spectra. 

Fig|4]shows the joint level spacing distribution for J = 
0, . . . , 4 and increasing number of central spins, where 
p{s) approaches closer and closer the Wigner surmise. 
To quantify this observation we use the ansatz 

Pis) = Bs^e-^'^ (2) 

with 

B = a(^((/? + 2)/a))'^+V(^((/3+l)/a))^+^ (3) 

A = {r{{l3 + 2)/a)/T{{P + l)/a)r (4) 

such that / p{s) — J sp{s) — 1. Clearly, a — 1, 
/3 — corresponds to an exponential distribution whereas 
a — 2, 13 = 1 reporduces the Wigner surmise. The 
above ansatz generalizes the Brody distribution given by 
a = P+1 [T31I3D]. As seen in Fig.[4j the numerical data is 
very well described by the above distribution, and for the 
case of 5 central spins the Wigner surmise of the GOE is 
almost reached. In particular, our level statistics do not 
show any odd/even-effects with respect to the number 
of central spins as predicted in Refs. [251 US]- We at- 
tribute this difference to the strongly restrictive assump- 
tions made there giving rise to additional symmetries. 

Moreover, in the case of two central spins it is instruc- 
tive to rewrite the Hamiltonian in the following form, 

n - {Si + s,)-Y,liA'i'^ + A<f^)n 

i 

+ iSi~S2)-Y.liA^^~A^^)n. (5) 

i 

The two central spins S — S1+S2 can couple to S* = 0, 1. 
Since the coupling to the singlet S — vanishes, the first 
line in Eq. ([5| is just the integrable Hamiltonian of a 
single central spin S = I, whereas the second line can 



be viewed as a perturbation. This term vanishes if the 
coupling constants are still random but chosen to be the 
same for each spin, A\^'' — \ resulting in an integrable 
model of two central spins, a prediction we have explicitly 
verified in our numerics; the latter model was also studied 
numerically in Ref. p5j . 

The models studied so far have a common spin bath, 
i.e. each bath spin couples without any further restric- 
tion to each central spin. Regarding the generic exam- 
ple of two neighboring quantum dot spin qubits this is 
not particularly realistic since in this geometry one can 
obviously identify groups of nuclear spins which couple 
strongly to one of the electron spins but weakly to the 
other. The extreme case is given by two separate spin 
baths where the central spins can be coupled via an ex- 
change interaction [311 132], 71' = Ti + JexSi ■ 82- Here 
we find numerically that even arbitrary small exchange 
parameters Jex break integrability and lead to level re- 
pulsion. The corresponding level spacing distributions, 
however, are less accurately described by the ansatz ([2|. 
On the other hand, for large \Jex\ the system approaches 
the integrable scenario since then the singlet and triplet 
subspace of the central spins are energetically more and 
more separated. 

Let us now discuss the influence of an external mag- 
netic field coupling to the central spins. In the case 
Nc = 1 the resulting model is known to be integrable 
[TOl [T7] , and also for A'c > 1 the Hamiltonian can still be 
represented as a real and symmetric matrix. Indeed, we 
have not seen any qualitative difference in the level spac- 
ing distribution with and without an external magnetic 
field. In particular, we have not found any sign for a tran- 
sition between the Gaussian orthogonal to the unitary 
ensemble (as appropriate for systems lacking time rever- 
sal symmetry [S]). In this sense, the application of an 
external magnetic field can be viewed as a "false symme- 
try breaking" which still preserves a "non-conventional 
time-reversal invariance" [5S1 [53] . We note that recent 
theoretical works predict different time dependencies of 
spin dynamics in different magnetic-field regimes |20| 121) . 
These observations are not reflected by the level statis- 
tics. Thus, decoherence and the occurrence of integra- 
bility or chaoticity are independent phenomena in such 
systems, at least as far as the role of magnetic fields is 
concerned. 

The data presented so far was obtained for bath spins 
of length / — 1/2. Motivated by the large nuclear spins 
in semiconductor materials, we have also performed sim- 
ulations for 1=1 which also do not show any qualitative 
difference to the previous case. This is indeed to be ex- 
pected since a spin bath of / = 1 can be obtained from 
a bath with J = 1/2 and twice the number of bath spins 
by grouping the spins into pairs and chosing the coupling 
parameters to be the same in each pair. Similar consid- 
erations apply to higher bath spins. 

Let us come back to the case of two central spins. As 
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seen in Figs. 2j4 this system appears to lie in between the 
integrable case and the predictions of random matrix the- 
ory. Thus, in the hght of the Bohigas-Giannoni-Schmitt 
conjecture [33], it is natural to speculate that the classi- 
cal counterpart of this system has a mixed phase space 
consisting of areas of regular and of chaotic dynamics. 
The classical limit of a quantum spin system is naturally 
obtained via spin-coherent states, and a pair of classical 
canonically conjugate variables p, q for each spin is given 
hy p = cos-d, q = ip, where ip are the usual angu- 
lar coordinates of the classical spin unit vector (35j . We 
have performed numerical Runge-Kutta simulations of 
such classical dynamics for Nc = 2 where one can easily 
treat systems of several thousand bath spins. However, to 
avoid the complications of such a high-dimensional phase 
space let us concentrate on the smallest nontrivial case of 
just two bath spins. Here we find indeed a close vicinity 
of regular and chaotic dynamics. An example is shown 
in Fig. [5] where we have plotted in the top panel a cut 



through the plane If 



as a function of p = S^, 



q = tan ^{Sf/Sf). The coupling constants are given by 
A^^' = 1 - Ai^^ = 0.75, ' = 1 - A[^'' = 0.52, and the 
initial condition is 5*1 = —5*2 = (0,0,-1), /i = (1,0,0), 
I2 = (0, 1,0). This arrangement leads obviously to very 
regular dynamics, in stark contrast with the bottom 
panel where we have used the same initial condition but 

introduced a minute change in one pair of coupling con- 

(2) (2) 

stants, A\ ' = 1 - A^^ = 0.5195, resulting in a clearly 
chaotic orbit with an inhoniogeneous phase space fill- 
ing. Note that the observation that certain phase space 
curves are overlaid in the figure is due to the fact that 
the remaining six phase space variables are not uniquely 
determined by the condition If = /| = and the con- 
served quantities 'H = 0,J=(1,1,0) but occur in several 
branches. 

In summary, we have investigated central spin mod- 
els via nearest-neighbor level spacing distributions. As 
the number of central spins increases a transition from 
Poissonian statistics to the Gaussian orthogonal ensem- 
ble sets in which can be described by a generalized Brody 
distribution. These observations are not affected by the 
finite system size in our numerical simulations and are 
unaltered upon applying an external magnetic field. In 
the transition region, the classical counterparts of the 
models studied have mixed phase space. 
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FIG. 1: (Color online) Nearest-neighbor level spacing for a 
system of a single central spin = 1/2 and 13 bath spins of 
length / = 1/2. The red curve is the exponential p{s) = e~^. 
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FIG. 2: (Color online) Nearest-neighbor level spacing for a 
system of a two central spins = S2 = 1/2 and 12 bath 
spins of length / = 1/2. The red curve is the GOE Wigner 
surmise. 




FIG. 3: (Color online) The same data as in the bottom right 
panel of Fig. [2] but for iV = 10 and = 11 bath spins. 




FIG. 4: (Color online) Joint level spacing distribution for 
J = 0, . . . , 4 and increasing number of central spins. The red 
dashed lines are the exponential function {Nc — 1) and the 
Wigner surmise {N^ > 1). The green solid lines are a fit to 
the generalized Brody distribution B. 




FIG. 5: A cut through the plane /f = /| = of two different 
phase space orbits demonstrating the close vicinity of regular 
and chaotic dynamics in a system of two central spins (see 
text). 



